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In this paper, we will generalize the Ganss-Bonnet gravity to an energy dependent Gauss-Bonnet 
theory of gravity, which we shall call as the Gauss-Bonnet gravity’s rainbow. We will also couple 
this theory to a Maxwell’s theory. We will analyze black hole solutions in this energy dependent 
Gauss-Bonnet gravity’s rainbow. We will calculate the modifications to the thermodynamics of 
black holes in the Gauss-Bonnet’s gravity’s rainbow. We will demonstrate that even though the 
thermodynamics of the black holes get modihed in the Gauss-Bonnet gravity’s rainbow, the hrst 
law of thermodynamics still holds for this modified thermodynamics. We will also comment on the 
thermal stability of the black hole solutions in this theory. 


I. INTRODUCTION 


The form of the standard energy-momentum dispersion relation is fixed by the Lorentz symmetry. Even though the 
Lorentz symmetry is one of the most important symmetries in nature, there are indications from various approaches 
to quantum gravity that the Lorentz symmetry might be violated in the ultraviolet limit Q-i- Thus, it is possible 
that the Lorentz symmetry is only an effective symmetry which holds in the infrared limit of quantum gravitational 
processes. As the standard energy-momentum dispersion relation depends on the Lorentz symmetry, it is expected 
that the standard energy-momentum dispersion relation will also get modified in the ultraviolet limit. In fact, it has 
been observed that such modification to the standard energy-momentum relation does occur in the discrete spacetime 
0, models based on string field theory Q, spacetime foam (jl, t he s p in- network in loop quantum gravity (LQG) 0, 
non-commutative geometry and Horava-Lifshitz gravitvpd| -m 

The modification of the standard energy-momentum dispersion relation has motivated the development of double 
special relativity [l^. In this theory, apart from the velocity of light being the maximum velocity attainable, there 
is also a maximum energy scale in nature. This energy scale is the Planck energy Ep, and it is not possible to 
for a particle to attain energies beyond this energy. The double special relativity has been generalized to curved 
spacetime, and this doubly general theory of relativity is called gravity’s rainbow [l4l| . In this theory, the geometry of 
spacetime depends on the energy of the test particle. So, particles of different energy see the geometry of spacetime 
differently. Hence, the geometry of spacetime is represented by a family of energy dependent metrics forming a 
rainbow of metrics. This is the reason the theory has been called gravity’s rainbow. In order to construct this theory, 
the modified energy-momentum dispersion relation is written as 


E^f\E/Ep)-p^g^{E/Ep) = m 


( 1 ) 


where Ep is the Planck energy. The functions f{E/Ep) and g{E/Ep) are called rainbow functions, and they are 
required to satisfy 


lim f(E/Ep) = 1, lim g{E/Ep) = 1. (2) 

E/Ep-^O E/Ep^O 

This condition is needed, as the theory is constrained to reproduce the standard dispersion relation in the infrared 
limit. Motivated by this energy dependent modification to the dispersion relation, the metric h{E) in gravity’s rainbow 
is written as [Tsil 

h{E) = g^\a{E)®eb{E). (3) 

Here the energy dependence of the frame fields can be written as 

^ fiEjEpf^' = giEjEpf^^ 
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where the tilde quantities refer to the energy independent frame fields. The energy at which the spacetime is probed is 
represented by E, and the maximum attainable energy is represented by Ep. So, if a text particle is used to probe the 
geometry of spacetime, then E is the energy of that test particle. Thus, by definition E cannot become greater than 
Ep [l5l| . As we will be applying gravity’s rainbow to the black hole thermodynamics, the energy E will correspond to 
the energy a quantum particle in the vicinity of the event horizon, which is emitted in the Hawking radiation 
It is possible to translate the uncertainty principle Ap > \/Ax into a bound on energy E > 1/Ax [^. In gravity’s 
rainbow, even though the metric depends on the energy, the usual uncertainty principle still holds [22j| . It has been 
demonstrated that the uncertainty in the position of the particle near the horizon should be equal to the radius of 
the event horizon radius 

A > 1/Ax « l/r+. (5) 

This energy E of the particle in the vicinity of the event horizon will be used in the rainbow functions. It may be noted 
that this energy A is a dynamical function of the radial coordinate [l^. We will not need the explicit dependence 
of this energy on the radial coordinate, but it is important to note that the rainbow functions are dynamical, and so 
they cannot be gauged away. Furthermore, this energy of a particle near the horizon is bounded by the Planck energy 
Ep, and it cannot increase to an arbitrary value. 

This bound on the energy also modifies the temperature of the black, and this modification to the temperature 
of the black hole can be used to calculate the corrections to the entropy of the black hole in gravity’s rainbow 0- 
The modification of the thermodynamics of black rings in gravity’s rainbow has also been calculated [2^ . It has been 
observed that the temperature of both the black holes and black rings starts to reduce after attaining a maximum 
value. Thus, at a critical size the temperature of the black holes and black rings becomes zero. At this value the 
entropy also becomes zero. Thus, the back hole stops radiating Hawking radiation, after reaching this critical radius. 
Thus, the gravity’s rainbow predicts the existence of a black remnant. It has been argued that such a remnant will 
form for all black objects [^. In fact, it has been explicitly demonstrated that such a remnant does form for Kerr 
black holes, Kerr-Newman black holes in de Sitter space, charged AdS black holes, higher dimensional Kerr-AdS black 
holes and black saturn [2^. 

It may be noted that the behavior of gravity’s rainbow depends on the form of the rainbow functions chosen. The 
construction of these rainbow functions is motivated from various theoretical and phenomenological considerations. 
Even though the form of these rainbow functions is different, the main property of these rainbow functions is that 
they make the geometry of spacetime dependent on the energy of the probe. This is the property which we will require 
in our analysis. So, in this paper, we will analyze the behavior of black hole solutions in the Gauss-Bonnet gravity’s 
rainbow using these phenomenologically motivated rainbow functions. 


II. GAUSS-BONNET GRAVITY’S RAINBOW 

It may be noted that it is possible to add higher order curvature terms to the original Lagrangian for general 
relativity. The theory obtained by adding quadratic powers of the curvature terms is called Gauss-Bonnet gravity. 
It is also possible to construct a Gauss-Bonnet-Maxwell gravity [2^-[2^. The Lagrangian of Gauss-Bonnet-Maxwell 
gravity can be written as 

Ltot = i? - 2A + aiRabcdR"’’"‘' - + R^) - T, (6) 

where Rated and A are the Riemann tensor and the cosmological constant, respectively. Here a denotes the Gauss- 
Bonnet (GB) coefficient, R = R'^ denotes the Ricci scalar and Red = Read denotes the Ricci tensor. The last term 
in Eq. © is the Maxwell invariant E = FabE°'^, where Fab = daAb — dbAa is the electromagnetic field tensor and 
At is the gauge potential. Now the fields equations for the the Gauss-Bonnet-Maxwell gravity Lagrangian m can be 
written as 

^ab + Agab + = --QabE + “^FacF^, ( 7 ) 

where 

VaE“'' = 0. 

Here the the Einstein tensor is denoted by Gf/ and 

G GB o / E? ryede o jy TDcd o jy jyc i jy jy \ 

ab — ^ [J^aede-^b ~ ^^aebd^ ~ -^^ac^b ' ^^ab) ~ 

i {RedefR^'^^^ - 4RedR‘^‘' + R^) gab. 


( 8 ) 
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In gravity’s rainbow the spacetime geometry depends on the energy of the probe {E). We can follow the methods 
used in the usual gravity’s rainbow [l3| , and absorb all the energy dependence of the Gauss-Bonnet gravity’s rainbow 
into the rainbow functions. So, we can write the energy dependent metric Gauss-Bonnet gravity’s rainbow as 


dr^ = —ds^ 




HE) 


1 



( 9 ) 


where 


d-2i-l 

dn^ =d9j+'^Y[ sin^ 9jdef (10) 

i=2 j=l 

Since we are looking for the black hole solutions with a radial electric field, we know that the nonzero components of 
the electromagnetic field are 


Ftr — —Frt- 


( 11 ) 


One can use Eq. m with the metric (|9l) to obtain the following explicit form of Ftr 

( 12 ) 

where q is an integration constant related to the electric charge of the black hole. In addition, it is easy to show that 
the following metric function satisfies all of the field equations o, simultaneously 


''W = ‘ + 2^ (‘- v’^) ■ 

with 

, ,, 8o' , (d-l)(ci-2)m (d-l)id-3)gy{E)^g{Ef 

9W = l+ (j_l)(d_2) 0+--- 

where m is an integration constant that is related to mass and a' = {d — 3)(d — 4 )q;. Taking into account Eqs. (fT^ 
and USD, it is clear that one cannot remove rainbow functions with rescaling. 

The choice of the rainbow functions f{E/Ep) and g{E/Ep) is very important for making predictions. That choice 
is preferred to be based on phenomenological motivations. Many proposals exist in the literature, we will some forms 
which have important phenomenological motivations. 

The rainbow functions motivated by the results obtained from loop quantum gravity and non-commutative geometry 
are given by [30j|-[3ll| 



f{E/Ep) = l, g{E/Ep) = ^l-r^{E/EpT. 
Now using these rainbow functions the metric is given by 


4'i(r) = 1 + 


(l - \/0i(’’)) : 


with 


0i(r) = 1 + 


8a' 


{d-l){d-2) 


A- 


( \ f 

f) 

r / \^ 

{d-l){d-2)m id-l){d-3)q'^\l-g(^^j 


2rd-i 


r,2d-4 


(15) 


(16) 


(17) 


The hard spectra from gamma-ray burster’s are given by can also be used to motivate the construction of rainbow 
functions [s^. These rainbow functions are given by 


PEjEp _ 1 

f{E/Ep)= , g{E/Ep) = l. 


( 18 ) 
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The metric corresponding to Eqs. (fT^ . (IT4t and (ITSl) . can be written as 


«'2(r) = 1 ^ - VQsW) . 


with 


02 (r) = 1 + 


8a' 


(d-l)(d-2) 


A- 


(d-l)(d-2)m (d-l)(d-3)q^E^(^e^ -ij 


V 


2j^d—l 


^2£^2j,2d-4 


(19) 


( 20 ) 


rainbow functions in which the velocity of light is constant have also been analyzed. This choice of rainbow functions 
is given by [3^ 


f{E/Ep)=giE/Ep) = 


1 


1-XE/Ep' 

The metric function corresponding to the case where the velocity of light is a constant can be written as m 

(i-fy 


( 21 ) 


^'3(?') = 1 + 


2a' 


(i - v'03(o) 


( 22 ) 


with 


03 (r) = 1 + 


8a' 


(d- l)(d-2) 


A- 


(d—l)(d—2 )to (d—l)(d —3)g^ 


V 


2rd—i 


r^2d-A 




(23) 


It may be noted that the bounds on the values of /3, rj, A has been analyzed using various theoretical and experimental 
considerations 1^. 


III. THERMODYNAMICS AND THERMAL STABILITY 


In this section, we will analyze the thermodynamics of the black hole solution obtained in the previous section. It 
is possible to obtain the Hawking temperature by using the surface gravity, 


T = 




2Ar^ 2f(E)(d-3)‘^q^ 
'¥{E) 


r+f{E) 


47r(d-2)(l + ^^^^) 


(24) 


According to the area law, the entropy of black holes is equal to one-quarter of the horizon area. This relation 
is valid for Einstein gravity, whereas we are not allowed to use it for higher derivative gravity. In this paper we 
investigate it in Gauss-Bonnet gravity’s rainbow. We can use the Wald formula for calculating the entropy 


S' = i / d^-^xy/^ (l -b 2aR^ 


(25) 


where R is the Ricci scalar for the induced metric jab on the (d — 2) dimensional boundary. We obtain 


-2 


S = 


Vd-2rX 

45'^-2(£;) 


1-b 


2(d-2)a'g^(E) 
(d — 4)r^ 


(26) 


which confirms that the obtained black hole solutions violate the area law. We should mention that Vd -2 denotes the 
volume of (d — 2)-dimensional sphere. 

Considering the Gauss law and calculating the flux of the electric field at infinity, one can find the electric charge 
of the black hole has the following form 


Q = 


Vd-2{d-3)q f{E) 
47r g‘^~^{E) 


(27) 
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We can calculate the electric potential $ as 

‘I’ = A^x^l 


- 


>r ->r+ „d-3' 


(28) 


where x^ is the null generator of the horizon. 

In order to calculate the finite mass of the black hole we will use the ADM (Arnowitt-Deser-Misner) approach for 
large values of r 


M = 


Vd -2 m{d-2) 
16^ f{E)gd-^E)- 


(29) 


Now using Eqs. (na), (iMi) and (12711 and considering M as a function of the extensive parameters S and Q, we obtain 
the following result, 


M{S,Q) = 


id-2) 


167rfiE)g‘i-^iE) 


g\E)ri-^ f 1 + ^X^ 


2Ar: 


d-1 


327r^ 


(d—l)(d—2) (d — 2)(d — 3)r'_f_ 

We can calculate the temperature and electric potential as the intensive parameter using the following relation 

dS 


(30) 


\dsjQ \dr+jQ/ \dr+jQ ’ \dQJs' 


(31) 


It may be noted that these relations are similar to the relations obtained in Eqs. (l24l) and (1281 Thus, the conserved 
and thermodynamic quantities satisfy the first law of thermodynamics. 


dM = TdS + 


(32) 


It may be noted that we can investigate the thermal stability of the charged black hole solutions of Gauss-Bonnet 
gravity’s rainbow through the canonical ensemble. In this ensemble, the electric charge is set as a fixed parameter 
and, therefore, the positivity of the heat capacity 






V j Q 


(33) 


is sufficient to ensure the local stability. Since we are investigating the physical black hole solutions with positive 
temperature, it is sufficient to examine the positivity of 


1 2 


Ti 


where 


fd'^M\ (rf-2)r+ [(d-4)ag2(^;) + _^ 

2gd-^iE)T2 


= -‘^-^E3lfS§pE + id-A){d-b)ag\E)Eg\E)rl 

(a 2jr_|_ 

2Ar^ 


(34) 


(d-2)(d-3)’ 


(35) 


To = — 


2(d - 3)q^PiE)g\E) Ud - 4)(2d - 7)agHE) + 


id-2) 


2d-8 


y \ 0 fi/ \ (d — 4)id — 9)ag^iE)r^ 

(d-4)2(d-5)aV(^^) + — \ j + 


2(d- 3) 


g^iE) [d-2 + 12aA (d - 4)] 
2{d-2)id-3) 


Krl 


(d-2)(d-3)2' 


(36) 
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FIG. 7: Model 1: Cq versus r+ for E = 1, Ep = 5, A = —1, q = 1, a = 2, n = 2, d = Q, and t] — 0 (solid line), t) — 2 (bold 
line) and rj = 5 (dotted line), ’’different scales” 





FIG. 8: Model 2: Cq versus r+ for E = 1, Ep = 5, A = —1, 5 = 1, a = 2, d = 6, and j3 = 0.1 (solid line), /? = 2 (bold line) 
and (3 = 5 (dotted line), ’’different scales” 


Thus, we are able to derive an explicit expression for the energy dependence of the heat capacity of a black hole in 
Gauss-Bonnet gravity’s rainbow. The positivity of this heat capacity ensures the local stability. 

In order to discuss the effects of gravity’s rainbow on thermal stability, we will analyze various different choices of 
rainbow functions. We plot heat capacity versus horizon radius for these different choices of rainbow functions. We 
also investigate the effects of d (spacetime dimensions), a (GB parameter) and the rainbow parameters on 

the heat capacity. 

The head capacity also depends on the dimensions of spacetime, as can be observed from Figs. [T][3] These figures 
demonstrate that divergent points appear in higher dimensions. Depending on the values of the parameters, one may 
find that for low dimensional cases, Cq is a smooth function and with a real positive root. However, by increasing 
the dimensions of spacetime, divergences can occur for Cq. In fact, two divergences occur and Cq is negative in 
between these two divergent points. This makes the system unstable. Although the global behavior of Cq is uniform 
for different choices of rainbow functions, their local properties such as slopes and the location of the root, along with 
the occurrence of divergences are different for different choices of rainbow functions. 

All figures indicate that Cq has a real positive root at r+o for all values of parameters. The location of r+o will be 
affected by variation of the parameters. It can be observed from Figs. fTUT^ that by increasing d, a and n the value of 
r_|_o decreases. Furthermore, by increasing 77 , (3 and A the value of r+o increases. It is possible to choose special values 
of the parameters such that the heat capacity does not diverge and remains a smooth function. However, it is also 
possible to fix the free parameters in such a way that Cq has two real positive roots at r+i and r _|_2 {r+i < 'r+ 2 )- In 
Figs. [TlfT51 we find that there is a critical value for the GB parameter (oc), such that a > Oc, and for this value of a, 
the heat capacity has two divergent points. Furthermore, Cq is a regular function for a < Uc- It may be noted that 
the value of Oc depends on the values of other parameters. In addition, considering the case of having two divergent 
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FIG. 12: Model 3: Cq versus r+ for E = 1, Ep = b, K = —1, q = 1, a = 2, d — G, and A = 0.1 (solid line), A = 1 (bold line) 
and A = 3 (dotted line), ’’different scales” 



FIG. 13: Model 1: Cq versus r+ for E = 1, Ep = 5, A = —1, g = 1, a = 4, t; = 2, d = 6, and n = 1 (solid line), n = 2 (bold 
line) and n = 3 (dotted line), ’’different scales” 


points in Fig. [TUI find that one can obtain a smooth Cq by increasing the value of 77 and (3 (decreasing A). The 
same behavior can be obtained by decreasing the value of n. In other words, since we regard 0 < -^ < 1, one finds 
that the effects of varying rj and (3 are opposite to the effects of varying n and A. Thus, the effect of gravity’s rainbow 
decreases by that decreasing 77 , f3 and increasing n, A. 


IV. CONCLUSION 

In this paper, motivated by the developments in gravity’s rainbow, we have developed a generalization of Gauss- 
Bonnet gravity to an energy dependent version of this theory. This new theory called Gauss-Bonnet gravity’s rainbow 
has been used for analyzing black hole solutions. In this theory, the metric in the Gauss-Bonnet gravity is made 
energy dependent in such a way that the regular Gauss-Bonnet gravity is recovered in the infrared limit. We have 
explicitly calculated the modifications to the thermodynamics of black holes in Gauss-Bonnet’s gravity’s rainbow. 
This was done by using various phenomenologically motivated rainbow functions. We used the rainbow functions 
whose construction was motivated by hard spectra of gamma-ray bursters at cosmological distances [s^. We also 
used the rainbow functions whose construction was motivated from results that have been obtained in loop quantum 
gravity and non-commutative geometry Finally, as it is possible to construct the rainbow functions in which the 
velocity of light is constant [s^, we also used these to analyze our results corresponding to such rainbow functions. We 
have also demonstrated that the first law of thermodynamics still holds for this modified thermodynamics. Finally, we 
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commented about the thermodynamical stability of these solutions. We showed that Cq has two different behaviors. 
In the first case it is a smooth and regular function with a real positive root at r+o • In this case black hole solutions are 
thermally stable for r+ > r+Q. For the second case, there are two singular (divergence) points for the heat capacity 
in which they are corresponding to second order phase transition points. We showed that variation of the model 
parameters, dimensionality and GB parameter can affect on the heat capacity behavior. 

It may be noted that in higher dimensions it is possible to construct various interesting solutions to black holes 
like black rings and black saturns. The black rings [2^, and black saturns [1^ have been analyzed using gravity’s 
rainbow. It has been demonstrated that in gravity’s rainbow these solutions also have a remnant. It will be interesting 
to construct black rings and black saturn solutions in Gauss-Bonnet gravity’s rainbow. It is expected that these 
solutions will also have a remnant because of the modification of the thermodynamics due to the rainbow functions. 
It has also been recently argued that the energy needed to create a mini black hole in a particle accelerator will 
increase in gravity’s rainbow. Thus, it would be possible to test a certain limit of gravity’s rainbow at the LHG 
[ 3 ^. This increase in energy occurs as the energy needed to form a mini black hole at the LHC has to be greater 
than the energy of the black hole remnant. It will be interesting to analyze a similar situation using Gauss-Bonnet 
gravity’s rainbow. It has been proposed that gravity’s rainbow can be used to address the black hole information 
paradox 36|-[33. So, it would also be interesting to repeat such an analysis for Gauss-Bonnet gravity’s rainbow. 
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